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ABSTRACT 

Line emission from dark matter is well motivated for some candidates e.g. sterile neutrinos. We 
present the first search for dark matter line emission in the 3 — 80 keV range in a pointed observation of 
the Bullet Cluster with NuSTAR. We do not detect any significant line emission and instead we derive 
upper limits (95% CL) on the flux, and interpret these constraints in the context of sterile neutrinos 
and more generic dark matter candidates. NuSTAR does not have the sensitivity to constrain the 
recently claimed line detection at 3.5 keV, but improves on the constraints for energies of 10 — 25 keV. 

Subject headings: Dark matter - line: identification — X-rays: galaxies: clusters 


1. INTRODUCTION 

Dark matter searches are a key pursuit of both astro¬ 
physics and particle physics. The scenario where dark 
matter is in a form of particles which provide gravity, but 
otherwise interact very weakly with ordinar y matter o r 
photons, is the most compelling (e.g. Taoso et aL| 2008). 
The most promising astrophysical objects tor searches 
for dark matter are clusters of galaxies, along with the 
Galactic Center and dwarf galaxy satellites to the Milky 
Way. Here, we consider the top end of the mass scale, 
galaxy clusters, with total masses - most of it in the form 
of dark matter - often exceeding 10 14 M 0 (~ 2 x 10 47 
g). Most clus ter mass estimates are inferred from X-ray 
observations (Vikhlinin et al. 2009 Mantz et al. 20101. 
Under the assumption of hydrostatic equilibrium. X-ray 
data provice a reasonably accurate model of the mass 
distribution, and imply cluster masses that are roughly 
consistent with masses measured via gravitational lens- 
ing, albeit both types of data are available for only a lim¬ 
ited number of objects. Their total masses are roughly 
five times greater t han the baryonic masse s infe rred from 


X-ray luminosities (von der Linden et al. 


2014). 


One possible particle candidate for dark m atter is 
the sterile neutrino (described in the reviews Boyarsky 
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. 20091 Kusenko 2009 Drewes|[2013[ and references 
link In the framework of this minimally neutrino 


et al 
therein) 

extension of the standard model (z/MSM), the lightest 
of the three sterile neutrinos provides the dark matter. 
The mass is basically unbound from theory, but some as¬ 
trophysical constraints apply. The mass is firmly bound 
from below through the phase space density of nearb y 


& Gunn 

1979 

) gives a model indej 

rendent 

of roughly U.4keV 

Boyarsky et al. 

2008b) 


can be increased if the production method is known; 
e.g., for reson ant production the bo undary is approxi¬ 
mately IkeV (Boyarsky et al. 2008b). An upper limit of 


a few hundred keV comes from a combination of produc¬ 
tion mechanisms and line emission searches. The ster¬ 
ile neutrino can, in principle, decay to two photons of 
equal energy, and thus the resulting signature would be 
a narrow emission line, corresponding to an energy of 
E 1 = m s c 2 / 2, where m s is the mass of the sterile neu¬ 
trino. The line width would be determined roughly by 
the velocity dispersion of dark matter particles, which 
for clusters is of the order of v ~ 1000 km sec -1 , which is 
smaller than the instrumental resolution of current ob¬ 
servatories. 

If we can limit the X-ray flux of a line at a specific en¬ 
ergy, and know the distance, we right away have a limit 
on the luminosity of the line from the cluster. Since we 
know the total mass of the cluster, we know how many 
sterile neutrinos there must be at a given assumed energy 
to provide the total mass of the cluster. Since we have a 
limit on the line luminosity, the ratio of the two is basi¬ 
cally the limit on the decay rate. A number of authors 
have reporte d searches for such emission lines in the soft 
X 


et al. 2006; 

Abazajian & Koushiappasl 

2006; Boyarsky 

et al. 2007b[ 

Riemer-Sprensen et al.12007 

Boyarsky et al. 


et "aT 


lew c 


I|20d9| |Loewenstein & Kusenko||2010[ |2012[) , with a 
laims of p oten tial detections that are yet to be con- 


yet 


firmed ([Lo cwenstein & Kusenko 2012 Boyarsky et al. 


2014 |Bulbul et al. 2014). f or an unambiguous detec- 


tion, the searches must avoid the spectral regions with 
line emission associated with atomic (or nuclear) transi¬ 
tions from the cluster gas, corresponding to any elements 
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w ith an expected appr eciable cosmic abundance. 


Bulbul et al. (2014) recently reported a possible sig¬ 


nature tor such a sterile neutrino at E 1 ~ 3.5 keV in 
stacked spectra of galaxy clusters observed with the 
XMM-Newton satellite; the result was confirmed at lower 


significance in a couple of individual clusters (Boyarsky 
et al. 2014), but it remains to be independently con- 
firrned in other types of dark matter objects, or using 
different instruments such as Suzaku jRiemer-Sprensen 
2014| Tamura et al.||2015[ ISekiya et al.||2015[ |Malyshev 
et al.|2014[|Anderson et ai.|2014|). Regardless, since there 
is no theoretical expectation as to the mass of the dark 
matter sterile neutrino beyond the broad rage described 
above, one should search for its signature in all accessible 
X-ray spectral bands. 

Here, we report on a search extending the energy range 
to the hard X-ray band, and thus the mass of the puta¬ 
tive sterile neutrino to twice the upper end of the NuS- 
TAR's bandpass, m s = 156 keV. Prior to the launch of 
NuSTAR , there were no sensitive spectral measurements 
beyond lOkeV, mainly b ecause sensitive meas urements 
require focussing optics (Harrison et al. 2013). While 
these energies have been searched tor li ne emission pre¬ 
viously using the cosmic background (Boyarsky et al. 
2006a 2008b|, this is the first search in a" pointed ob¬ 


servation with focussing optics. 

Recently the NuSTAR team observed and reported the 
results on one well-studied cluster - the “Bullet Cluster” 


(Wik et al. 2014). There are many previous observa- 
tions covering its X-ray and lensing properties (Marke- 
vitch et ah] 20021 12004| [Clowe et al. 2006[ |Paraticz et al. 
2012|). This galaxy cluster, at z = 0.296, is perhaps 
best known for the detailed comparison of the distribu¬ 
tion of dark matter as inferred from gravitational lensing 
to the X-ray emitting gas. The lack of spatial align¬ 
ment betwe en the two distributions reported in |Clowe | 
et al. ([2006j ) on t he basis of a weak lensing analysis and 


Bradac et al. (2006) using a joint weak and strong 


lensing analysis, is often considered to be one of the 
strong est arguments for the existence of dark matter par¬ 
ticles. Bradac et al. (2006) estimate a total mass of the 
cluster to be 5 x 10 i4 M 0 within the central 500 kpc. Due 
to the offset between the mass and X-ray emitting gas, 
the Bullet Cluster provides an excellent low-background 
environment for dark matter searches. 

The Chandra X-ray data of the Bullet have been 
searched for isolated X-ray emission lines out to 10 keV 


dRiemer-Sprensen et al.|2007| |Boyarsky et al.|2008c|) but 
here we extend the ran ge fd~C 80keV by using the NuS- 


Wik et al. 


TAR observations from 
Our approach (described in bee. 13.11) is to fit the data 
with an adequate model describing the emission by the 
hot gas in the cluster, and then search for the improve¬ 
ment to the fit by adding an isolated emission line of 
varying energy but with fixed width (determined by the 
instrumental resolution). The detection (or limit on the 
flux of the line) provide the measurement (or limit) of the 
flux and thus luminosity of decay photons, yielding the 
volume-integrated decay rate of the putative dark matter 
particle (presented in Sec. [4]). 

2. OBSERVATIONS 

The Bull et Cluster was o bserved by NuSTAR in two 
epochs (see Wik et al.pOM for details) for a total com¬ 


bined exposure of 266 ks. To filter the events, stan¬ 
dard pipeline processing (HEASoft v6.13 and NuS- 
TARDAS vl.1.1) was applied along with strict criteria 
regarding passages through the South Atlantic Anomaly 
(SAA) and a “tentacle”-like region of higher activity near 
part of the SAA; i.e. in the call to the general processing 
routine that creates Level 2 data products, nupipeline, 
the following flags were included: SAAM0DE=STRICT and 
TENTACLE=yes. No strong fluctuations are present in 
light curves culled from the cleaned events, suggesting 
a stable background, so no further time periods were ex¬ 
cluded. 

From the cleaned event files, spectra and response files 
were created using nuproducts. The call to nuproducts 
included extended=yes, most appropriate for extended 
sources, which weights the response files based on the dis¬ 
tribution of events within the extraction region, assuming 
that to be equivalent to the true extent of the source. 
We extracted spectra for each of the regions shown in 
Fig. 0 (details in Table IT]). The regions were chosen to 
maximise the amount ofaark matter within the field of 
view while minimising gas emission. The “Peak” region 
contains the leading mass peak but excludes the shock 
front. It is id entic al to the analysed region in |Boyarsky| 
et al.| (|2008c|) and |Riemer-Sprensen et al.| (|2007|), while 
for the “Half Peak” region, the extracted area is larger 
to compensate for the lower spatial resolution of NuS¬ 
TAR compared to Chandra. The “Left” region contains 
the trailing mass peak with the exclusion of the main gas 
emission. For each region we chose an offset background 
region of identical shape at a location with similar gas 
emission but less mass as inferred from the lensing map. 
The resulting spectra for all three regions are shown in 

Fig- [2} 

3. ANALYSIS 
3.1. Spectral modelling 

We fit the spectra of both observations and both de¬ 
tectors simul taneo usly using the Xspec spectral fitting 
package (Arnaud 1996) and explored two different ap¬ 
proaches to background modelling and subtraction. The 
spectra were binned to contain at least three counts per 
bin, and we assume the background to be Poisson dis¬ 
tributed and consequently us e Cash-stat istics ( C ) to op¬ 
timise the parameter values (Cash 1979). 

In the first method of background treatment we sim¬ 
ply subtract the spectrum of the dark matter offset region 
from the spectrum of the source region. Since the spectra 
are almost identical we can fit any residuals with a single 
power law (the statistics for each region are given in Ta¬ 
bled]). Subsequently, we added a Gaussian to represent a 
single emission line at a fixed energy and flux to the best 
fit model above, and searched for the improvement of C 
as a function of line energy and intensity, allowing for 
simultaneous variation of the power law parameters and 
the Gaussian normalisation. We consider line intensities 
from 0 to 10 -5 photons cm -2 sec -1 and line energies of 3 
to 80 keV in steps of A E = 0.1 keV. We assume that the 
line is narrow compared to the detector resolution, fixing 
the intrinsic line width at 0.001 keV, and noting that as 
long as the assumed width is less than ~ 0.03 keV, our 
results do not change. 

The reduction of the C parameter by the extra line is 
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Figure 1 . Image of the NuSTAR observations of the Bullet Cluster, with the X-ray contours from Chandra overlaid in black (|Riemer-| 
S 0 rensen et al.|2007|), and the weak lensing contours in green (|Clowe et al.|2006|). The white circles illustrate the source regions with the 
second circle excluded to avoid gas emission. The magenta circles are the background regions chosen to have similar gas emission as the 
source regions, but much less mass. 


Table 1 

The first section provides the coordinates of the regions illustrated in Fig. [T] and the dark matter mass within each field of view based on 
weak gravitational lensing. The second section contains the fit statistics for each of the regions for the methods of background subtraction 
and modelling. The spectra were binned to a minimum of three counts per bin before the analysis. 


Region 

Peak 

Half Peak 

Left 

Peak+Left 

Included center (RA, Dec) [degrees] 

104.56825, -55.941758 

104.56825, -55.941758 

104.64978, -55.951826 

— 

Included radius [aremin] 

0.66 

0.66 

1.00 

— 

Excluded center (RA, Dec) [degrees] 

104.58827, -55.942086 

104.59071, -55.994209 

104.62326, -55.944488 

— 

Excluded radius [aremin] 

0.375 

0.66 

1.00 

— 

Mass [Mq] 

5.70 x 10 13 

4.46 x 10 13 

1.02 x 10 14 

1.57 x 10 14 

Power law c-parameter/dof 

1460.3/1355 

1195.6/1105 

1967.4/1889 

3405.2/3221 

Background model c-parameter/dof 

1610.0/1646 

1411.6/1347 

2380.3/2315 

3990.3/3933 



Figure 2. The spectra of the three different regions: Peak (blue), 
Half Peak (green), Left (red). The individual exposures have been 
stacked and rebinned for visualisation purposes. The bump around 
30keV comes from instrumental lines and only affects the back¬ 
ground modelling method, not the background subtraction. 

shown as a function of line energy in the upper panel of 
Fig. [3] At most energies, the additional Gaussian does 
not lead to any significant improvement (|AC| < 9 for 
1105-3821 degrees of freedom), and instead we constrain 
the flux by increasing the Gaussian normalization and 
refitting all other parameters until AC = C — Cb aS e = 
2.71, corresponding to the one-sided 95% confidence level 
marginalised over the power law normalizations. These 
flux levels are shown in Fig. [3] for the Peak and Left 


regions as well as for the combined analysis. 

In the second approach we mo del the back ground in¬ 
stead of subtracting it. In Wik et al. (20141 the back¬ 
ground emission in the Bullet/ Cluster was thoroughly 
investigated and we use their results as a baseline model 
for the background, and check if there is room for any line 
emission above this model. The model consists of four 
components: i) the aperture background (smooth gradi¬ 
ent across the detector with a normalisation uncertainty 
of 10%); ii) a focused cosmic ray background from unre¬ 
solved sources (smooth with a normalisation uncertainty 
of 10%); and iii) instrumental continuum (we use a 10% 
normalisation uncertainty even though the systematic 
uncertainty is probably much smaller); iv) the thermal 
solar continuum and instrumental lines from reflections 
(smooth component with a normalisation uncertainty of 
10% plus known detector emission li nes). Additionally 
we fit a line free plasma model (apec, Smith et al.|[20 01) 
to account for any gas emission from the cluster. The 
redshift and abundances of the plasma model are kept 
fixed at z = 0.296 and A = 0.2 respect ively (consistent 
with the value found in (Wik et al. 2014)), while the tem¬ 
perature is required to be the same tor both detectors and 
both observations, but with individual normalisations. 

The 3—150 keV energy interval is well fitted by the 
background model with the statistics given in Table [l] 
As before, for each energy between 3 and 80 keV in steps 
of A Ej = 0.1 keV we add a Gaussian and determine the 
best fit normalisation considering line intensities from 
—10 -5 to 10 -5 photons cm -2 sec -1 . The line intensities 
are allowed to be negative to account for overestimation 
of the background. The normalisation is then increased 
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Figure 3. Upper panel: The C statistics improvement (negative) from adding an extra Gaussian at the corresponding rest frame energy 
(Ery = (1 + z)E ohs ). Lower panel: The individual derived flux limits (95% CL) for the Peak region (blue) as well as the Peak and Left 
regions combined (thick black). The solid lines show the results of modelling the background of the nearby dark matter offset region, 
and the dotted lines show the results of subtracting the nearby dark matter offset region. The results from the individual regions are all 
consistent with each other and the constraints tighten by combining the regions. The background modelling generally provides slightly 
stronger constraints, but with gaps where the preferred flux is negative. 


until AC 2 = 2.71. We derive upper flux limits from the 
Gaussian alone as well as from the base model plus the 
Gaussian integrated over the full width at half maximum 
(FWHM) of the NuSTA R spect ral resolution crudely ap¬ 
proximated by ( Harrison et al.||2013 l 

( 1 ) 


A Ei 


FWHM 


= 0.01.EL + 0.3 keV. 


The addition of the Gaussian only improves the model 
by |AG| > 9 around 50 and 84keV where AG ss —12, 
reflected in weaker constraints at those energies. While 
|A|C > 9 would appear to imply a significant detection, 
we need to take the look-elsewhere effect into account 


(Gross & Vitells 2010). The energy of the line is un- 
known and by scanning over energy we perform a num¬ 
ber of independent searches which increases the chance of 
seeing statistical outliers. Consequently the probability 
of detection is degraded by the number of independent 
attempts (given the spectra resolution of NuSTAR we 
search of the order of 150 independent energies). 

Line-like features may arise from fluorescent and ac¬ 
tivation induced instrumental lines if imperfectly mod¬ 
elled, or due to statistical/systematic fluctuations be- 
tween ta rge t and background regions (see appendix of 
Wik et al. 2014). These lines are strongest between 
20 — 30keV, and may explain the region of larger C- 


values at those energies. 

The flux limits from the two approaches for the back¬ 
ground subtraction are compared in Fig. [3j The solid 
lines show the results of modelling the background of the 
nearby dark matter offset region, and the dotted lines 
show the results of subtracting the nearby dark matter 
offset region. The results from the individual regions are 
consistent with each other, and the constraints tighten 
by combining the regions. The background modelling 
generally provides slightly stronger constraints, but with 
gaps where the preferred flux is negative. The results are 
very similar, and we consider the background treatment 
to be robust. We also derive flux limits from fitting the 
two regions simultaneously, taking their mass difference 
and consequently expected signal strength into account. 

While the background modelling provide slightly 
stronger constraints, there is a risk of spurious line de¬ 
tections from small residuals due to a gain drift between 
the actual observation and the time of collection of data 
from which the background model was constructed. This 
is unlikely to be the case for direct background subtrac¬ 
tion. 


3.2. Mass within field of view 
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The weak lensing shear mapi^from Clowe et al. (2006) 
can be integrated to provide the masses within the three 
regions for a fiducial cosmology of O m = 0.3, Da = 0.7, 
Ho = 70 km sec -1 Mpc -1 . The map contours are shown 
in Fi g. [1] and the obtained region masses are given in 

Table ET 


Paraficz et al. (20121 presented a mass map of the Bul¬ 
let Cluster based on strong lensing rather than weak lens¬ 
ing. This map provides region masses that are almost 
twice as big as for the weak lensing map. This indicates 
an uncertainty on the mass estimates of the order of 50%. 
In the remaining analysis an d plots we conservat ively use 
the values from the map by Clowe et al. 2006]) and as¬ 
sume that the entire mass is made up of dark matter 
since the observed regions have been chosen to minimise 
the gas presence, so for the Peak region M gas ss 0.1M to t- 

4. RESULTS 

4.1. Sterile neutrinos 

Assuming all of the dark matter to be sterile neutri¬ 
nos, we can interpret the flux constraints from Fig. [3] 
in terms of the sterile neutrino parameters of mass, m s 
and mixing angle, sin 2 (20), where the latter describes the 
mixing probability with the lightest of th e active neutri¬ 


nos. The constraints are converted as (|Riemer-S0rensen 
et al. |[2006j Boyarsky et al.|[2007b ): 


sin 2 (20) < 


10 


18 


/ Fobs \ 

/ rn s \ 

r(M fov /M 0 )i 

^ erg cm -2 sec -1 J 

V keV/ 

(Da/M pc) 2 _ 


( 2 ) 


Again, the NuSTAR constraints do not significantly 
improve on existing constraints, but being a pointed ob¬ 
servation with robust background treatment they provide 
an important cross-check on previous analysis of diffuse 
emission observed with HEA P, INTEGRAL, COMP- 


TEL and EGR ET (Yiiksel & Kistler 2008 Boyarsky et al. 
2008bl|2006a). 


5. DISCUSSION 

5.1. Possible detections 

Unfortunately, the claims of possible line detections 
(Loewenstein & Kusenkol 2010 Bulbul et al. 2014 Bo- 
yarsky et ai.|2014|) all lie below the lower sensitivity cut- 
ott tor AluSTAR of 3.89 keV introduced by the redshift of 


the Bullet Cluster. 


5.2. Possible improvements 

The Milky Way halo has provided stro ng constraints 
on line emission in the X-ray ran ge (|Riemer-S0rensen 
et al. 2006 Riemer-Sprensen 2014). The advantage of 
local constraints is that the dark matter source is nearby 
and there exists a wealth of observations, but the disad¬ 
vantage is the number of sources of non-thermal “back¬ 
ground” radiation and the uncertainty of the inner mass 
profile. The background radiation issue can be reduced 
significantly by point source removal if one has sufficient 
spatial resolution. This is now becoming possible with 
NuSTAR , and will be investigated further in future work. 
The profile uncertainty problem can be mitigated some¬ 
what by excluding the centre of the halo from the anal¬ 
ysis. 


where F 0 b s is the observed flux limit, Mf ov is the total 
dark matter mass within the field of view, and Dl is the 
luminosity distance. 

As illustrated in Fig. [4j very large mixing angles will 
lead to over-production of dark matter and are conse¬ 
quently ruled out. Similarly, the resonant produc tion 


mechanisms require a primordial lepton asymmetry (Bo¬ 


yarsky et al. 2008b), that may affect Big Bang nucleosyn¬ 
thesis and cosmic element abundances. This gives a lower 
limit on the mixing angle. As mentioned in Sec. [l] the 
mass ra nge is limited from below by the TremaineAjunn 
bound (|Tremaine & Gunn||1979| |Boyarsky et al.||2008b[ 
and from above by line emission searches. The NuSTA 
Bullet Cluster constraints from the Peak region alone are 
weaker than previously existing constraints from the dif¬ 
fuse cosmic background, but by combining the Peak and 
Left regions we get similar constraints. This provides an 
important independent cross-check of the many assump¬ 
tions about, e.g., source distribution that goes into the 
diffuse background constraints. 


4.2. Generic dark matter constraints 

In Fig. [5] we present the constraints on generic dark 
matter decays leading to photon emission for two photons 
per decay. For one-photon interactions, the constraints 
are weaker by a factor of two. For particles of Majorana 
type where the particles are their own anti-particles, the 
interaction probability doubles and the constraints are 
strengthened by a factor of two. 

12 http://flamingos.astro.uf1.edu/le0657/index.html 


6. SUMMARY 

We have searched NuSTAR observations of the Bullet 
Cluster for exotic line emission over the 3 — 80 keV in¬ 
terval. No significant line flux was found and we have 
derived upper limits on the possible line emission flux. 
While the constraints are similar to previous constraints 
from the cosmic background emission, this is the first 
time a search has been performed in this energy interval 
using a pointed observation. The constraints can be im¬ 
proved by longer observations or usin g diff erent target s, 
e.g., the Milky Way halo (as in Riemer-Sprensen]2014). 


This research made use of data from the NuSTAR mis¬ 
sion, a project led by the California Institute of Tech¬ 
nology, managed by the Jet Propulsion Laboratory, and 
funded by NASA, and it also made use of the NuSTAR 
Data Analysis Software (NuSTARDAS) jointly devel¬ 
oped by the ASI Science Data Center (ASDC, Italy) and 
the California Institute of Technology (USA). Facilities: 
NuSTAR. 
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Tremaine-Gunn 
■ Thermal production 

Maximal lepton asymmetry 
M31 Chandra 

— Chandra Sgr A* 95% CL, NFW 
HEAO diffuse background 
INTEGRAL diffuse background 
NuSTAR Peak+Left subtraction, 95% CL 
NuSTAR Peak+Left modelling, 95% CL 
Boyarsky et al. 2014 
O O Bulbul et al. 2014 


m s [keV] 


Figure 4. The derived constraints on sterile neutrino mixing angle as a function of mass compared to a selection of previous constr a ints. 
The upper grey region is excluded by over-production of dark matter, and the lower by Bie Bane nucleosynthesis (jL ainc & Shaposhnikov| 
|2008||Canetti et al.|2013|. Below m s = 1 keV the Tremaine-G unn l imit applies ( Tremaine & Gunn 1979 Boyarsky ct al. 20U8bj. me green 
points indicate the pot ential signal from i Bulbul et al.l2014ll and iBovarskv et al -liitl 111 an ' ' ' 

Chandra (Horiuchi et al.l2 014). Chandra Sgr A* |Riemer-borcnscn 2014|) 
liffuse background 'JTioyarsky et al.|2008b[|. The Nuijl'Ati LSuiiet Cluster 


pom 

bounds: M31 
INTEGRAL diffuse 
the diffuse cosmic background, but provide an important cross-check. 


and the coloured shaded regions arc obser vation 
""HEAO diffuse background i|Boyarsky et al. |200Ga[l , 
constraints are similar to previous constraints trom 



E 1 [GeV] 


INTEGRAL, COMPTEL, EGRET 
INTEGRAL diffuse background 
HEAO diffuse background 
Chandra Sgr A* 95% CL 
M31 Chandra 

Prototype cryogenic spectrometer 
Chandra LETG NGC3227 
NuSTAR Peak+Left subtraction, 95% CL 
NuSTAR Peak+Left modelling, 95% CL 


Figure 5. The derived constraints on two-photon decay rat es as a function of photon energy compared to previous constraints: IN¬ 
TEGRAL, COMPTEL and EGRE T (lYuksel & Kistlcr 20081, INTEGRAL d iffuse background (Boyarsky et al. 2008b), HEAO diffuse 
background (Boyarsky ct al. 200 6all. Chandra Sgr A * ijfuemcr-Sdrcnscn 2014), M31 Cha ndra (Horiuchi ct al. 11 4(1. prototype cryogenic 
spectrometer ( [Boyarsky et al. |2DP7a) , Chandra LETG NGG3227 (|Bazzocchi et al.|2008|). These constraints applies to all Majorana type 
dark matter candidates with mono-energetic photon emission in the relevant energy range. The NuSTAR Bullet Cluster constraints are 
similar to previous constraints from the diffuse cosmic background, but provide an important cross-check based on different assumptions. 































































































NuSTAR line emission constraints 


7 


Boyarsky, A., Iakubovskyi, D., Ruchayskiy, O., & Savchenko, V. 
2008a, MNRAS, 387, 1361 

Boyarsky, A., Malyshev, D., Neronov, A., & Ruchayskiy, O. 
2008b, MNRAS, 387, 1345 

Boyarsky, A., Neronov, A., Ruchayskiy, O., & Shaposhnikov, M. 
2006a, MNRAS, 370, 213 

Boyarsky, A., Neronov, A., Ruchayskiy, O., Shaposhnikov, M., & 
Tkachev, I. 2006b, Physical Review Letters, 97, 261302 
Boyarsky, A., Nevalainen, J., & Ruchayskiy, O. 2007b, A&A, 471, 
51 

Boyarsky, A., Ruchayskiy, O., Iakubovskyi, D., & Franse, J. 

2014, Physical Review Letters, 113, 251301 
Boyarsky, A., Ruchayskiy, O., & Markevitch, M. 2008c, ApJ, 673, 
752 

Boyarsky, A., Ruchayskiy, O., & Shaposhnikov, M. 2009, Annual 
Review of Nuclear and Particle Science, 59, 191 
Bradac, M., Clowe, D., Gonzalez, A. H., et al. 2006, ApJ, 652, 937 
Bulbul, E., Markevitch, M., Foster, A., et al. 2014, ApJ, 789, 13 
Canetti, L., Drewes, M., & Shaposhnikov, M. 2013, Physical 
Review Letters, 110, 061801 
Cash, W. 1979, ApJ, 228, 939 

Clowe, D., Bradac, M., Gonzalez, A. H., et al. 2006, ApJ, 648, 
L109 

Drewes, M. 2013, International Journal of Modern Physics E, 22, 
30019 

Gross, E., & Vitells, O. 2010, European Physical Journal C, 70, 
525 

Harrison, F. A., Craig, W. W., Christensen, F. E., et al. 2013, 
ApJ, 770, 103 

Horiuchi, S., Humphrey, P. J., Onorbe, J., et al. 2014, 

Phys. Rev. D, 89, 025017 
Kusenko, A. 2009, Phys. Rep., 481. 1 

Laine, M., & Shaposhnikov, M. 2008, Journal of Cosmology and 
Astro-Particle Physics, 6, 31 
Loewenstein, M., & Kusenko, A. 2010, ApJ, 714, 652 
—. 2012, ApJ, 751, 82 


Loewenstein, M., Kusenko, A., & Biermann, P. L. 2009, ApJ, 700, 
426 

Malyshev, D., Neronov, A., & Eckert, D. 2014, Phys. Rev. D, 90, 
103506 

Mantz, A., Allen, S. W., Rapetti, D., & Ebeling, H. 2010, 
MNRAS, 406, 1759 

Markevitch, M., Gonzalez, A. H., Clowe, D., et al. 2004, ApJ, 

606, 819 

Markevitch, M., Gonzalez, A. H., David, L., et al. 2002, ApJ, 567, 
L27 

Paraficz, D., Kneib, J.-P., Richard, J., et al. 2012, ArXiv e-prints, 
arXiv: 1209.0384 

Riemer-Sdrensen, S. 2014, ArXiv e-prints, arXiv:1405.7943 

Riemer-Sdrensen, S., & Hansen, S. H. 2009, A&A, 500, L37 

Riemer-Sdrensen, S., Hansen, S. H., & Pedersen, K. 2006, ApJ, 
644, L33 

Riemer-Sdrensen, S., Pedersen, K., Hansen, S. H., & Dahle, H. 
2007, Phys. Rev. D, 76, 043524 

Sekiya, N., Yamasaki, N. Y., & Mitsuda, K. 2015, ArXiv e-prints, 
arXiv: 1504.02826 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, 

J. C. 2001, ApJ, 556, L91 

Tamura, T., Iizuka, R., Maeda, Y., Mitsuda, K., & Yamasaki, 

N. Y. 2015, PASJ, 67, 23 

Taoso, M., Bertone, G., & Masiero, A. 2008, Journal of 
Cosmology and Astroparticle Physics, 3, 22 

Tremaine, S., & Gunn, J. E. 1979, Physical Review Letters, 42, 
407 

Viklilinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009, ApJ, 
692, 1060 

von der Linden, A., Allen, M. T., Applegate, D. E., et al. 2014, 
MNRAS, 439, 2 

Wik, D. R., Hornstrup, A., Molendi, S., et al. 2014, ApJ, 792, 48 

Yiiksel, H., & Kistler, M. D. 2008, Phys. Rev. D, 78, 023502 



